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ABSTRACT 

We study the emission and dissipation of acoustic waves from cool dense clouds in pressure equilibrium 
with a hot, volume-filling dilute gas component. In our model, the clouds are exposed to a source of ionizing 
radiation whose flux level varies with time, forcing the clouds to pulsate. We estimate the rate at which acoustic 
energy is radiated away by an ensemble of clouds and the rate at which it is absorbed by, and dissipated in, the 
hot dilute phase. We show that acoustic energy can be a substantial heating source of the hot gas phase when 
the mass in the cool component is a substantial fraction of the total gas mass. We investigate the applicability 
of our results to the multiphase media of several astrophysical systems, including quasar outflows and cooling 
flows. We find that acoustic heating can have a substantial effect on the thermal properties of the hot phase in 
those systems. 

Subject headings: cooling flows — ISM: clouds — ISM: jets and outflows — galaxies: Seyfert — quasars: 
general — waves 
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1. INTRODUCTION 

Multi-thermal-phase gas configurations are often encoun- 
tered in astrophysics. For example, interstellar gas is known 
to consist of several thermally distinct components whose 
temperatures are in the range 10 2 - 10 6 K, (e.g, McKee & 
Ostriker 1977). Gaseous halos of L* galaxies constitute cool 
(10 4 K) gas condensations embedded in a hot (10 6 K) volume- 
filling virialized gas (e.g., Mo & Miralda-Escude 1996, Che- 
louche et al. 2007 and references therein). The intracluster 
medium (ICM) of nearby clusters shows rich structure of cool 
ionized filaments embedded in a hot X-ray emitting material 
and extending to large scales (e.g., Conselice et al. 2001, 
Fabian et al. 2003). A similar structure is also observed 
around radio-loud quasars (RLQs) where a considerable frac- 
tion of the mass of the gaseous nebulae lies in cool clumps of 
gas (e.g., Crawford & Fabian 1992, Fu & Stockton 2006). 
More recently, it has been realized that many astrophysi- 
cal outflows are also multiphase with cool condensed mate- 
rial embedded in a hotter and more dilute ambient medium. 
These systems include: planetary nebulae (e.g., Meaburn et 
al. 1992), stellar winds (e.g., Bouret et al. 2005), and quasar 
flows (e.g., Das et al. 2005). 

The physical mechanisms responsible for multiphase gas 
configurations in different systems are quite diverse. A par- 
tial list includes thermal instability and compression by radia- 
tive shocks in hot material, and evaporation and gravitational 
binding of cool material. While our understanding of the com- 
plex physics leading to multiphase structures is incomplete, it 
is well established that such systems are a wide-spread phe- 
nomenon. 

Recent studies suggest that some multiphase systems seem 
to suffer from a heat "deficit". One example is that of highly 
ionized X-ray outflows from the central regions of Seyfert 
galaxies that are believed to be thermally driven and are seen 
to reach high velocities despite the effect of adiabatic cooling 
(e.g., Chelouche & Netzer 2005). An analogous case to that 
are the large scale narrow line region (NLR) outflows in active 
galactic nuclei (AGN) where the effect of adiabatic cooling is 
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even more pronounced (e.g., Everett & Murray 2007). These 
discrepancies may be alleviated if a yet unidentified heating 
source is lurking in those systems which would balance adi- 
abatic losses. Another example is that of cooling flows in 
which the time for radiative cooling of the hot component ap- 
pears to be short compared to the Hubble time, yet the gas 
does not cool significantly. Various explanations have been 
put forward to explain the cooling flow problem including 
thermal conduction (e.g., Zakamska & Narayan 2003), PdV 
work done by expanding radio bubbles that are inflated by (re- 
cently activated) AGN (e.g., Fabian et al. 2003 and references 
therein; see also Mathews et al. 2006), and the dissipation 
of gravitational energy via cool condensations falling through 
the hot gas (Murray & Lin 2004; see also Dekel & Birnboim 
2007). 

In this paper we wish to address the general problem of the 
heating/cooling balance in multiphase media around sources 
with time-varying ionizing fluxes. Flux variability is a char- 
acteristic of many astronomical objects and of active (accret- 
ing) systems in particular [see e.g., Geha et al. (2003) for 
the case of quasars and Zezas et al. (2007) for a study of X- 
ray binaries]. While the examples given in this paper focus 
on quasar-related phenomena, we emphasize that the general 
processes outlined here may apply to other systems such as 
X-ray binaries and variable stars. 

This paper is organized as follows: In §2 we present a gen- 
eral analytic formalism that describes the radiative forcing of 
clouds which causes them to pulsate and emit acoustic waves 
into their surrounding medium. Absorption of acoustic ra- 
diation and several heat transport mechanisms are also dis- 
cussed. We address a few specific systems to which this phys- 
ical mechanism may apply in §3. Numerical calculations con- 
firming some aspects of our analytic approximations are also 
presented. Our conclusions follow in §4. 

2. THE EMISSION AND DISSIPATION OF ACOUSTIC WAVES 

2.1. General Setup 

Here we consider a generic configuration of a multiphase 
medium by assuming dense and cool gas condensations 
(clouds) with time-averaged temperature T c embedded in a 
hot, dilute, volume-filling medium at temperature 7), (see Fig. 
1). We assume that the two phases are, in a time-average 
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FIG. 1. — A schematic (not-to-scale) view of the multiphase medium con- 
sidered in this work. Clouds (dark circles) are embedded in a hot dilute 
medium and are roughly in pressure equilibrium with it (see text). The clouds 
are assumed to be spherical with a radius R c , density, p c , and temperature, T c 
(see inset). The entire volume is illuminated by an ionizing source of varying 
flux (upper-right corner) that radiatively forces the clouds to pulsate and emit 
acoustic waves into their surrounding medium. 

sense, in pressure equilibrium such that 

T 

Ph = Pc^r, (1) 

where p c , Ph are the densities of the cool and hot compo- 
nent, respectively. Characterizing such an equilibrium re- 
quires knowledge of the relevant heating and cooling mech- 
anisms in each phase. These processes are complicated, but 
we note that, generally, cooling is °< p 2 ; hence the cooling 
time would be shorter for the dense and cool phase compared 
to the dilute and hot one. In a steady-state thermal balance, 
this is also true for the heating rates. 

Below we show that a fraction of the photon energy which 
is absorbed by the cool gas can be redirected to heat the hot 
ambient medium by means of acoustic waves. Such a pro- 
cess naturally occurs when the cool phase is photo-heated 
by a time-varying radiation field (e.g., that of a quasar). 
Clouds distributed throughout the volume of the system would 
change their pressure due to the varying radiative heating and 
would therefore pulsate. As such, they would act as bells (or 
monopole loudspeakers) emitting acoustic radiation into their 
environment. Unlike Murray & Lin (2004), our model for 
the heating of the volume-filling hot medium does not require 
differential velocities between the phases to operate, and does 
not necessarily lead to the grinding of the cool phase. 

2.2. Emission 

Consider an isolated spherical cloud which is exposed to a 
source of ionizing flux that varies with time t (Fig. 1). Con- 
sider also a specific mode of oscillation such that the flux, 

F(t)=F Q + 8Fe" ot , (2) 

where f ffl = 2n(0~ l and 8F is the amplitude of variation. We 
further assume that 

max(4 0l ,4 at ) < t m < max( 

^cool ! Aieat ) (3) 

with C>1> f heat ( f cooi> O bein 8 the cooling time and heating 
time of the cool (hot) phase, respectively (see §3 for the dis- 
cussion of specific systems). In this regime, the cloud reacts 
instantaneously to the varying flux by changing its tempera- 
ture, and hence its pressure (we check the validity of this con- 
dition for specific systems in §3). So, for example, as the level 



of ionizing flux rises, more energetic electrons are injected 
into the gas thereby increasing its temperature. The cloud, be- 
ing over-pressured compared to its environment, would then 
expand into the ambient medium. For large enough flux vari- 
ations (8F/F ~ 1), the cloud's temperature and pressure (P) 
variations satisfy AT C /T C ~ AP/P ~ 1 (see §3 for detailed pho- 
toionization calculations). In this case, the cloud will expand 
at roughly its sound speed, c c . We neglect any time depen- 
dence of c c in the present analysis as this is proportional to 
y/ 'AP/P which is of order unity. We note that since t ra <C f<? ool , 
the hot phase will not react to flux variations of this duration 
and its properties may be considered constant over t m . The 
last condition may be relaxed if the hot component is in a 
range of parameter space where <9ln(r) /d\\\(F) — 0, as would 
be the case for gas at the Compton temperature. 

In what follows, we shall work in the long wavelength 
regime, namely, t m > R c /ch where R c is the cloud radius (see 
below for a more general definition of t a which takes into ac- 
count the finite cooling time of the cloud). At much higher 
frequencies a phase develops between the gas pressure and 
velocity (as they become orthogonal on small scales) so that 
the cycle-averaged radiation is small. 

Acoustic flux is often over-estimated in the literature. When 
a cloud — or a piston — expands subsonically (indeed c c <C c/, 
by construction in our case), the ambient medium in its im- 
mediate environment would be effectively incompressible and 
would pulsate with the cloud not producing sound (i.e., com- 
pression) waves. Sound waves are produced once the sur- 
rounding medium has not had time to respond to the changing 
volume of the region interior to it which occurs at a distance 
~ Chta from the cloud. This distance marks the transition to 
the radiative zone (see appendix). From the continuity condi- 
tion, the velocity at that location 

<*> 

Hence, for large enough f ffl , | -C c c . The acoustic luminosity 
radiated away by a single pulsating cloud would therefore be 

^acoustic = 2%R 2 c p h C 2 c C h ^1 + ^j^j • (5) 

Clearly, for large enough t a , L acoU stic becomes very small. 

It is interesting to compare the total energy emitted over 
one cycle to the energy stored in the hot gas within a shell 
defined by the surface of the cloud at a distance r = R c and 
r = R c + Cht a . As | declines rapidly with r, we can write 

^acoustic ^G) Rc 

4/lnR? c p h c 2 c " c~^,' 

Hence, for long enough wavelengths, this ratio is <C 1 and 
most of the PdV work done by the cloud by expanding into 
the ambient medium is returned to it during the compression 
phase (e.g., Leighton 1997 and the appendix). This is also 
why typical loudspeakers are rather inefficient in producing 
sound converting, on average, only of the electric power 
to audible acoustic power. 

The ratio of acoustic energy emitted by all the clouds over 
the cooling timescale of the hot phase (from adiabatic or ra- 
diative losses) and the thermal energy stored in the hot phase, 
Eh, is 

^acoustic _ 3 £y / \ 2 c /i f cool (y c /i f (o \ 2 ,~ 
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with By being the volume filling factor of the cool phase. The 
last expression is valid in the limit Cht a /R c 3> 1 and e„ -C 1. 
Clearly, a pre-requisite for effective acoustic heating is that 
^acoustic > Eh- This is achieved for short t m . Nevertheless, t m 
cannot be arbitrarily short since it must satisfy 
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r ffl >^max(^. ^ , 
\c h A(T,p c ,F / 



(8) 



where A is the net cooling/heating rate per unit volume aver- 
aged over half a cycle. If the converse is true then the clouds' 
thermal state would change little with time and radiative forc- 
ing would be less efficient (this is shown to be the case in §3). 
For t w > t%, Elastic/Eh oc f^ 2 , hence t% is the mode which is 
most efficient in converting photon flux to acoustic flux. 

It is interesting to compare the acoustic luminosity emitted 
by the cloud to its photo-heating rate 



^acoustic ^ Rc 

An/'iRlA ~ c h t m 
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At a 
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(9) 



Clearly, only a small fraction of the photo-heating luminosity 
is transformed to acoustic luminosity (e.g., at most ~1% will 
be converted to acoustic luminosity for T c = 10 4 K and 7), = 
10 6 K). 

For the case in which the emitted acoustic radiation is fully 
absorbed in the system then acoustic heating would overcome 
the cooling of the hot phase once 



e v A(p c ) ( c c 



l-e v A(p h ) 
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(10) 



where we have used equation 9 and our assumption of pres- 
sure equilibrium. If the ratio of the cooling rates for the hot 
and cold phases is ~ (p/,/ p e ) 2 (as would be the case if, for ex- 
ample, the thermal properties of the hot phase are determined 
primarily by photo-heating; see §3) then we require that the 
mass fraction of the cool component, Em, satisfies 



Pc 

£ M ~ — £v 
Ph 



1. 



(11) 



This brings us to an important conclusion: the mass of the 
cool phase should be at least comparable to that of the hot 
phase for acoustic heating to be effective. We discuss this 
result more quantitatively in §3.1-3.3. 

It is possible to analytically estimate the typical size of 
clouds that can give rise to significant acoustic heating. A nat- 
ural scale comes from the requirement for a minimal t % which 
occurs for R c /ch = p c c^ /A (see equation 8 and the numerical 
results of §3). This gives the following order-of-magnitude 
estimate 



R, 



cool ^ 

l/c h ' 



(12) 



i.e., R c /l is of order the ratio of the cooling time of the cool 
phase to the sound-crossing timescale of the hot phase. A 
specific case worth mentioning is that of adiabatically cool- 
ing/heating, thermally-driven systems (such as X-ray and 
NLR quasar outflows). Here, the thermal state of the hot 
medium satisfies 

P~A— . (13) 

Ch 



That is, the time it take for the entire volume to heat (or cool) 
is of order the dynamical time. Using this property of such 
systems, we can readily estimate estimate R c to be, 



Rr 



I 



I 



(14) 



[in this case, the required number of clouds in the system is 
roughly (p c / p/,) 5 ]. We emphasize, however, that this provides 
only a rough estimate for R c and that, in fact, there is a range 
of cloud sizes which can induce significant acoustic heating, 
as we demonstrate for a few specific systems in §3. 

2.3. Absorption 

The emission of acoustic energy by a cool cloud does not 
guarantee that the energy will dissipate over the relevant spa- 
tial scales required to effectively heat up the volume. Dissipa- 
tion of acoustic energy in the hot medium occurs via absorp- 
tion (and conductance) as well as by steepening of acoustic 
waves into weak shocks. The latter process is negligible in our 
case since the clouds pulsate subsonically and the wave am- 
plitude is quickly diminished in the spherical case (c.f., Stein 
& Schwartz 1972). 

The relevant dissipation length scale due to absorption of 
radiation is (e.g., Mathews et al. 2006) 



h — M~ 



(15) 



where fi is the effective viscosity (including both the ef- 
fects of viscosity and conduction; e.g., Fabian et al. 2005). 
Here we take the usual Braginskii-Spitzer (e.g., Spitzer 
1962) value for high temperature solar composition gas ~ 



5 II 

\0~ ib r\T h gcm^'s -1 (e.g., Lang 1999) which is applica- 
ble to e.g., the ICM for 77 = 10~ 2 - 1 (Narayan & Medvedev 
2001, Fabian et al. 2005; but see Ettori & Fabian 2000 for 
possibly much lower values). Due to the lack of observational 
constraints on the conductivity of other systems (e.g., quasar 
outflows), we shall assume values similar to those of the ICM. 

Rewriting equation 7 with dissipation included, we obtain a 
generalized condition for the importance of acoustic heating, 



^dissipated 
acoustic 



(l-e-'/fc) 



>1. 



(16) 



For Id ^> /, absorption is inefficient and most of the acoustic 
energy (a fraction <~ 1 — (I /Id) of it) will escape the system. 
Absorption is efficient over the entire volume for lj ~ /. For 
Id -C /, absorption is very efficient yet is localized to the en- 
vironment of the clouds. In this case several additional con- 
ditions are required to effectively heat up the entire volume. 
These are discussed below. 

2.3.1. The case of Id <Si I 

For volumetric heating to be important, we require many 
clouds which are homogeneously distributed throughout the 
volume of the system and/or the presence of some heat trans- 
port mechanism. For clouds distributed randomly in space 
(as indeed seems to be observationally supported; see §1), the 



characteristic distance between clouds is 



,-1/3 



R c . If no 



heat transport mechanism is operating then we require 

1/3 h_ . 
Ev R r ~ l 



(17) 
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so that dissipation of acoustic waves due to all clouds occurs 
throughout the volume. A much less restrictive case is that in 
which heat conduction and/or advection are present. In this 
case, even if the above condition is not satisfied, heat could 
still be distributed over the entire volume. 

Heat conduction is due to free electrons whose velocity is, 
v e ~ y / 3k B T e /m e ~ 6 x 10 3 -\/7/ ! /10 6 K km s _1 . The propa- 
gation of electrons through a plasma takes the form of a ran- 
dom walk with some mean free path, X e . The largest uncer- 
tainty arises in estimating X e , which depends on the strength 
and configuration of the magnetic field, B, as well as on the 
pressure and density of the (hot) medium. If B = then 
X e = Xj? =0 ~ 10 4 r /l 2 /(p/ ! /m//) cm (where m# is the mass of 
a hydrogen atom; see Cowie & McKee 1977). When B > 
the medium becomes anisotropic with X e being of order the 
Larmor radius across magnetic field lines. For a detailed dis- 
cussion of these issues see Lazarian (2006, and references 
therein). Here we assume that 10~ 2 < X e /X^ =0 < 1 (Narayan 
& Medvedev 2001). Hot electrons must traverse a distance of 
order the typical cloud separation over ?p 00 i, so that heat could 
be uniformly distributed throughout the volume. This results 
in the following condition for efficient conduction: 



p -2/3 ff 2 



< 1. 



(18) 



In addition to conductance, heat may be advected by e.g., 
turbulent motion in the gas (we do not consider here the ad- 
ditional heating by the dissipation of turbulent energy; e.g., 
Chelouche & Netzer 2005). For a Kolmogorov-type transonic 
turbulence, the velocity of eddies of size l e is v e = Ch(l e /l) 1 ^ 
and for efficient advection we therefore require 



c /i f cool 

2.4. Caveats 



(19) 



We have outlined a general analytic framework for study- 
ing the emission of acoustic waves by pulsating cool, dense 
clouds and the dissipation of such waves in a hot surround- 
ing medium. We have not addressed the issue of cloud 
longevity, which could be important as clouds are prone to 
various hydrodynamic instabilities such as Kelvin-Helmholtz 
and Rayleigh-Taylor, as well as evaporation/condensation. 
Our understanding of these processes is rather limited in the 
general astrophysical context as well as in the specific systems 
considered below, and we make no attempt to further address 
this here. That said, the ubiquity of multiphase structures in 
astrophysical systems implies that such configurations are, at 
least statistically, quasi-steady-state phenomena. 

A further caveat concerns the role of magnetic fields in de- 
termining the propagation of sound waves in the medium. 
Better estimates for the magnitude of (magneto-)acoustic 
heating will require not only detailed numerical simulations 
but also a better understanding of the individual systems to 
which this process may apply. 

3. APPLICATIONS TO ASTROPHYSICAL SYSTEMS 

In this section we examine how acoustic heating can be im- 
portant in astrophysical environments. It is not the purpose of 
this paper to cover all possible cases but rather to study the ap- 
plicability of the process to a few specific systems. In all cases 
considered below there is either an active ionizing source with 




FIG. 2. — Time-dependent calculations of the temperature of a cloud which 
is exposed to a time varying flux of a sinusoidal form (denoted by a dot- 
ted line; see legend). All quantities are shown relative to their mean val- 
ues. t a was varied in the range (0.1 — 100) x /£ ool as calculated from the 
steady-state solution. As expected, the cloud thermal state follows the flux 
variations for t a 2> f£ j . A phase develops and the temperature variation am- 
plitude is somewhat lower for t a ~ (f j . The cloud's temperature variation is 
quenched for t a <g t£ ool . These calculations justify our approximation of the 
heating/cooling rates by their steady-state values (see §3). 

a varying flux or one has recently been active. Here we con- 
sider examples pertaining to quasar outflows and to cooling 
flows. To this end, we have taken fiducial parameter values 
from the literature concerning the physical properties of the 
cool and hot phase but note the appreciable uncertainties as- 
sociated with them. These parameter values include rough 
estimates for T c , Tf,, and p c . Estimates for the total size of 
the system, I, and its cooling time, f* j are also taken from 
the literature. The characteristic values of ey and R c are very 
poorly constrained in those systems and are therefore treated 
as free parameters. 

A major difficulty when estimating the response of the 
cloud to radiative forcing is to properly assess the net 
heating/cooling function A (equation 8) which depends on 
whether the gas is in photoionization and thermal equilibrium 
at all times (as would be the case if the heating and cooling 
timescales are the shortest in the problem) or whether its ther- 
mal state depends on its history in some intricate way. To keep 
the problem tractable, we estimate the net heating/cooling 
function by the value of the cooling function under steady- 
state conditions (in which case its values match those of the 
heating function). For solar metalicity, optically-thin gas ex- 
posed to a type-I quasar continuum (Sazonov et al. 2005), 
the heating rate does not have a strong dependence on the 
temperature, varying by a factor of order unity in the range 
10 4 < T c < 10 6 K, so that an adequate approximation for our 
purpose is 



A(r,p) ~ 10 2i (p c /m H ) 2 erg s cm" 



(20) 



This value shall be used to obtain an order of magnitude esti- 
mate for the heating/cooling timescale of the cloud (equation 
8). 

To check the validity of the above approximation for A 
in cases where strict thermal equilibrium is not justified, we 
have calculated several time-dependent photoionization mod- 
els using CLOUDY c07.02 (Ferland et al 1998) and assum- 
ing a light-curve, of a fixed amplitude, varying the mode 
frequency, i.e., ta>- Ischoric conditions are assumed for the 
calculations since c t /c/, <C 1 and we focus on the thresh- 
old of the long wavelength regime (see below). Calculations 
were carried out for the case of p c /niH = 10 5 cm -3 and plac- 
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(b) Extended NLR Outflows 

FIG. 3. — Contour level plots of the ratio of acoustic heating to cooling as a function of e v and R c (normalized to the size of the system, /) for the systems 
considered in this paper (§3.2-3.4). For any R c , it is assumed that there exists a mode with t a = t% for which acoustic heating is evaluated. Contours with values 
> 1 indicate regions in the parameter space where acoustic heating dominates over cooling in the hot phase. Three astrophysical systems are considered: X-ray 
outflows in Seyfert I galaxies (left), extended NLR outflows in Seyfert galaxies {middle), and cooling flows (right). Color coding of contours is the same for all 
panels. We note that the contours' shape near £y = 1 is somewhat uncertain due to the relatively little volume occupied by the hot material and the importance 
of absorption and scattering of acoustic waves from the cool phase in this case. Clearly, for all systems there exist a wide range in the parameter space where 
acoustic heating is important (or even dominant; see text). This occurs when the cool gas mass is comparable to or larger than that of the hot component. Smaller 
clouds are more effective in heating the system as long as R c /ch is greater than their cooling time. For quasar outflows, the results for the cloud size and the 
minimal £y, where acoustic heating becomes important, are consistent with analytic scalings (see §3) and some observational constraints (see text). 



ing the cloud at 1 pc from the quasar so that its steady-state 
temperature — if it were exposed to a constant flux F — is 
i~ 8 x 10 4 K (corresponding roughly to X-ray outflows tem- 
peratures; see §3.1). The results, shown in Figure. 2, indicate 
that our estimate for the net cooling function (i.e., ffLj) us- 
ing steady-state values is indeed reasonable. In particular, the 
cloud does not react to the varying flux level of the source 
for f co <C 

^cooi anc ^ ra diative forcing would be inefficient. The 
cloud temperature does follow the light-curve for t m ~ ^ ool , 
albeit with a phase, and its temperature fluctuations, AT C , are 
almost as large as those obtained for t m 3> £? oo] (where little 
or no phase is present since the cloud reacts instantly to the 
source's flux level). We note here that AT C depends not only 
on AF but also on dln(T) / dln(F): in particular, for the cho- 
sen ionizing continuum, large flux variations are expected for 
T c ~ a few x 10 4 — 10 5 K. The proper treatment of marginally 
stable or thermally unstable regions (Krolik & Rriss 2001) is 
beyond the scope of this paper. 

In the above, we have not addressed the possibility of 
changes to the spectrum of the ionizing source which are 
known to occur (e.g., Netzer et al. 2003) and could influ- 
ence the gas temperature even when the total integrated flux 
remains unchanged. Such a case can be easily accommodated 
within our general formalism (§2) by considering F to be the 
relevant ionizing flux level (e.g., that just above the Lyman 
edge for T c ~ 10 4 K gas or around soft X-ray energies for 
T c ~ 10 5 K gas). As spectral variations of active sources are 
not well constrained by observations and their inclusion will 
add little to our understanding of the general effect, we choose 
to not address this issue here. 

Quasars are known to vary considerably {8F/F ~ 1) over a 
wide range of timescales, from minutes to years (e.g., Giveon 
et al. 1999, Netzer et al. 2002 and references therein). While 
current observations do not directly constrain the power of 
such variations on the timescales and wavebands most rele- 
vant to our study, it is reasonable to assume that the ionizing 
flux of quasars varies considerably over the entire relevant 
frequency range. This motivates us to adopt the following 
approach: we assume that for any choice of system parame- 
ters (e.g., cloud size, density, and temperature), there exists 
a mode in the quasar light-curve oscillating with t a ~ and 
having a large amplitude of flux variation {8F/F ~ 1) . 



Our calculations indicate that, for all systems considered 
here, uniformly distributed heating by acoustic waves is pos- 
sible so that one (or more) of equations 17 — 19 is always 
satisfied when the more restrictive, equation 16 is. Generally, 
smaller clouds (shorter t m ) occupying a large fraction of the 
volume give rise to more efficient heating. Nevertheless, as 
clouds cannot be infinitely small, on account of their cooling 
time being longer than t m (equation 8), there is a interme- 
diate range of cloud sizes with a total volume above a cer- 
tain threshold, which is most efficient in heating the ambient 
medium. The details depend on the system in question as we 
now discuss. 

3.1. X-ray Outflows in Seyfert 1 Galaxies 

The parameters for X-ray gaseous outflows from Seyfert 
1 galaxies which we adopt here are: T c = 8 X 10 4 K, 7), = 
10 7 K, p c /m H = 3 x 10 5 cm" 3 , / = 3pc and ^ ! qol = 10 4 years 
(Chelouche & Netzer 2005). Considering equation 16, there is 
a large range of the parameter space in which acoustic heating 
is considerable and can balance adiabatic cooling (see Fig. 
3). In particular, such heating would be effective - that is, 
comparable to radiative heating - for £y > 10~ 2 , i.e., once the 
mass in the cool phase is comparable or larger than that of 
the hot phase (c.f. equation 11). Cloud sizes should be in 
the range 10 14 — 10 15 cm (R c /l ~ 10~ 4 ; see equation 14) for 
acoustic heating to be effective. Such values are consistent 
with full occultation of the X-ray source in Seyfert 1 galaxies 
(e.g., Chelouche & Netzer 2005). We note that if £y ~ 0.1 
(i.e., Em > 10) then acoustic heating may even dominate over 
photo-heating resulting in a hot-phase temperature that would 
exceed the Compton temperature. 

3.2. Extended NLR Outflows in Seyfert Galaxies 

Little is known about the physical properties of NLR flows, 
and the existence of a hot medium is mainly inferred from 
the dynamics of the low temperature, optically-detected gas. 
We take fiducial values of T c = 10 4 K, T h = 10 6 K and con- 
sider multiphase gas on ~ 30 pc scales (e.g., Das et al. 2005, 
Everett & Murray 2007). For luminosities typical of type- 
I Seyfert galaxies, the implied density of the hot phase on 
such scales is of order 10cm -3 (Everett & Murray 2007; see 
also Chelouche & Netzer 2005) hence p c /niH = 10 3 cm~ 3 . 
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Here, the dynamical time is of order the adiabatic-cooling 
time which is of order 10 s years (Everett & Murray 2007). 
In this case we find a wide range in parameter space where 
acoustic heating can be important. Everett & Murray (2007) 
quote values for the filling factor of the cool phase as high 
as 10~ 2 . Figure 3 shows that, in this case, the emission of 
acoustic energy by clouds with R c ~ 10 16 cm could contribute 
significantly to the heating-cooling balance of the hot phase. 
Clouds of this size cannot be individually resolved by cur- 
rent observations yet similar cloud sizes have been inferred 
from density and column density measurements of UV ab- 
sorber properties on similar scales (e.g., Rraemer et al. 2001 
and references therein). Arguments applying to X-ray out- 
flows (equation 18,19) apply here too. 

3.3. Cooling Flows in Galaxy Clusters 

Observations of RLQs show the presence of intrinsically 
bright O II and O III emission line regions with apparent sizes 
of a fewxlOkpc that consist of numerous unresolved cool 
clouds (e.g., Fu & Stockton 2006 and references therein). 
It is thought that at least some RLQs and radio galaxies in- 
habit the cores of galaxy clusters that possess massive cool- 
ing flows toward their centers (e.g., Crawford & Fabian 1992). 
In this case, the optical emission is thought to originate from 
cool clouds embedded within the hot medium (Crawford & 
Fabian 1992). It is not known whether such multiphase gas 
configurations are a property of all cooling flows in galaxy 
clusters; nor is it clear if all cooling flow clusters which do 
not cool effectively contain active or recently active nuclei 
(though they do seem to have a > 50% duty-cycle for jet ac- 
tivity; Fabian & Sanders 2006). Nevertheless, we can use 
typical cloud parameters inferred from optical observations 
and ambient medium parameters deduced from X-ray studies: 
T c = 2x 10 4 K, T h = 10 7 K, p c /m H = 10 cm" 3 , l = 5x 10 4 pc 
and f^ ! ool = 7 x 10 9 years (see e.g., Dunn & Fabian 2006 and 
Fu & Stockton 2006). Figure 3 shows that for e v > 10~ 3 (in- 
dicating comparable mass in the hot and cool phase; c.f. Fu & 
Stockton 2006), clouds of size ~ 0.01 — 1 pc could efficiently 
balance radiative cooling in those systems. Interestingly, sim- 
ilar cloud sizes have been deduced by Hamann et al. (2001) 
for the clouds seen in absorption toward quasars. Larger val- 
ues of £y (hence Sm) would be required for efficient acous- 
tic heating if the combined duty-cycle of intracluster active 
sources is much smaller than unity (Shen et al. 2007). Given 



the uncertainties associated with the duty-cycle of quasars in 
cooling flow environments and the properties of the cooling 
flow itself (e.g., concerning its multiphase structure), it is dif- 
ficult to assess just how important acoustic heating is, in the 
form considered here. 

4. CONCLUSIONS 

In this paper we have presented a simple model for estimat- 
ing the effect of acoustic heating by radiatively-forced pulsat- 
ing clouds on the thermal state of a hot ambient medium. We 
find that smaller clouds are more efficient acoustic emitters as 
long as their cooling time is longer than the oscillation period. 
Acoustic heating is significant when the mass of cool clouds 
is comparable to, or larger than that of the hot phase. In par- 
ticular, for large enough mass in cool gas, acoustic heating 
may even dominate over photo-heating. This implies the pos- 
sible existence of higher temperature gas around sources with 
variable flux compared to non-variable objects with similar 
characteristics. This work demonstrates that heating of mul- 
tiphase photoionized gas may depend not only on the mean 
flux of the source but also on deviations from it, as well as 
on the structure of the medium itself. We show that there ex- 
ists a range of clouds sizes which is most efficient in emitting 
acoustic energy. Applying this model to quasar outflows and 
cooling flows, we find that acoustic heating could be impor- 
tant in these cases if the typical size of cool clouds is of order 
10~ 4 the size of the system. Better understanding of the mul- 
tiphase nature of photoionized environments in terms of the 
relevant densities, cloud sizes, and cooling rates, as well as 
the properties of the ionizing source (i.e., its power spectrum 
of variations) are required to properly assess this effect in the 
broader astrophysical context. 
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APPENDIX 

Lagrangian perturbations for the equation of motion of the hot component in spherical coordinates gives, 

a / 1 d 



#S=c£V(V.S) 



c h 



-A 



(1) 



whose general solution is 

t J = [ah{kr) + by l {kr)]e im , (2) 

where j\, y\ are the spherical Bessel functions. The boundary conditions for the problem are: i) a forced oscillator at the wall 
of the cloud satisfying E, (r = R c ,t) = c t exp(;'of ) and ii) a radiatively expanding wave solution at infinity such that E, (r — > °°) oc 
exp(;©£ — kr). A solution is obtained upon matching the boundary conditions and is shown in Figure 4. Clearly, the flux amplitude 
drops quickly, and by a large factor, from the surface of the cloud (on the left side) toward the radiative zone. Hence, most of the 
energy emitted by the pulsating cloud at the start of the cycle is returned to it toward the end so that only a small fraction of the 
energy is acoustically radiated away. The acoustic flux at large distances from the cloud surface is constant in this example since 
no dissipative mechanism was included in equation Al. 



